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I. INTRODUCTION 


The aerodynamics of long rod finned projectiles has been surveyed 
by many investigatorsl>2,3,4,5,6,7,8 with principal efforts directed to 
resolving the individual contributions of the component hardware (nose, 
body and tail) in the determination of the static aerodynamic coefficients. 
Only rarely has the trajectory problem been addressed, although the 
specific analysis of the dispersion of the projectile at the target 
has been examined carefully in the open literature^*^®. Previous searches 
for predictive algorithms considering the complete free-flight package 
were not rewarding; and this led to the issue of several additional 


^C.H. Murphy 3 "Free Flight Motion of Symmetric Missiles"^ BEL Report No. 
1216^ July 1963 (AD #442757). 

9 

AMCP 706-280j "Design of Aerodyncmioally Stabilized Free Rockets", 1968. 

3 

W.F. Braun, "Aerodynamic Data for Small Arms Projectiles", BRL Report 
No. 1630, January 1973 (AD #9097571). 

i 

4 * 

H.W. Liepmann and A. Roshko, Element of Gasdynamics, John Wiley and 
Sons, Inc., New York, 1957. 

^A.H. Shapiro, The Dynamics and Thermodynamics of Compressible Fluid 
Flow, Volume I, The Ronald Press Company, New York, 1953. 

g 

M. Piddington, "The Aerodynamic Characteristics of a SPIW Projectile", 
BRL Memorandum Report No. 1594, September 1964 (AD #355679.).. 

7 

E.R. Dickenson, "Some Aerodynamic Effects of Blunting a Projectile 
Nose", BRL Memorandum Report No. 1596, September 196.4 (AD #451977)... 

8 

L.C. MacAllister, "Drag and Stability Properties of the XM144 Flechette 
with Various Read Shapes", BRL Memorandum Report No. 1981, May 1969 
(AD #854724). 


Q 

W.J. Gallagher, "Elements Which Have Contributed to Dispersion in the 
90/40 mm Projectile", BRL Report No, 10.13, March 1957 (AD 3135.3061, 

10 

J.D. Nicolaides, C.W. Ingram, "Analysis of the Jump and Dispersion of 
Flechettes", Prepared for U.S. Army, Frankford Arsenal under Contract 
No. DAAA 25-71-C0447. 
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constituant reports » > » » which, however, had the ultimate 

objective of composing a more comprehensive publication which would 
enhance estimating facility for this particular class of flight projectile. 

This report is an attempt to scribe the envelope of the preliminary 
design configuration. The drag, static moment and normal force coefficients 
are determined from the simplified geometry of the projectile and the 
bounding aerodynamics conditions; i.e., sea level flat-fire, 2<M<5 and 
low yaw as developed in the previous reports^^’^^ with the accuracy factor 
(related to the dispersion) and retardation similarly extracted 
The program is written in Tektronix 4051 Basic language with graphic 
extensions. The program flowchart, sample graphics screen display and 
program listing are given in Appendices A, B, and C, respectively. 

The example provided for illustration has not been extensively 
flight tested; indeed, a complete set of data to include range aerodynamics, 
retardation and dispersion for a given projectile over the selected 
Mach excursion is simply not available. The influence of the dynamic 
aerodynamic coefficients, the damping and the Magnus, on the flight 
performance is even more difficult to certify. However, what disparate 
data that is published does suggest that the given criteria for design 
are viable. It is also obvious that local modifications to the program, 
based on individual experience, are certainly inevitable and are 
actively encouraged. 


11 

W.F. Donovarij "Procedure for Estimating Zero Yaw Drag Coefficient for 
Long Rod Projectiles at Mach Numbers from 2 to 5", BEL Memorandum 
Report No. ARBRL-MR-02819, March 1978 (AD U0S4326). 

12 

W.F. Donovan^ "One Factor Affecting the Dispersion of Long Rod 
Penetrators"^ BRL Memorandvm Report No. ARBRL-MR~02846^ June 1977 
(AD #A058696). 

13 

W.F. Donovan, "Simplified Determination of Retardation for Kinetic 
Energy Projectiles", BRL Memorandum Report No. ARBRL-MR-02994, 

February 1980 (AD #083299). 

^^W.F. Donovan, "Algorithm for Estimating Aerodynamic Static Moments 
of Long Rod Penetrators at 2<M<5, " BRL Memorc^ion Report No. ARBRL-MR- 
03020, May 1980 (AD #086095). 

^^W.F. Donovan, "Hypothetical Zero Yaw Drag Coefficient of Kinetic 
Energy Projectiles Between M = 5 and M = 10", BRL Memorandum Report 
No. ARBRL-MR-03041, August 1980 (AD # 090009). 


10 



II. PROCEDURE 


The analytic expressions for the various parameters are taken from 
the previous reports and tabulated. Recent test results have been used 
to modify the equations, where indicated, to improve correspondence 
between test and prediction. The fin-body interference factor has been 
adjusted to recognize a fin masking effect of the boundary layer and 
a separate drag contribution has been assigned to the presence of the 
driving grooves. Caliber notation, where a representative linear 
dimension provides a reference length and the mass/force dimension is 
converted to a specific gravity equivalent, is employed extensively. 

This rather novel practice allows a direct comparison of projectile 
performance on the basis of a "ballistic weight" quite independent of 
physical dimension. 

Figure 1 is an outline of a typical long rod projectile. Tables 
1, 2 and 3 contain the equations used in developing the programs. The 
initialization instructions are presented in Appendix A and the full 
program listings follow. The program provides an outline plot of the 
projectile with the physical properties and gives a complete chart 
description of the aerodynamic coefficients with accuracy factor, 
decrementing velocity and initial yaw cycle estimates. The accompanying 
print-out lists the discrete component contributions and the retardation 
over the expected range of the flight to M = 5. 

III. RESULTS 

The 25/12 projectile, recently tested in the BRL aerodynamics range, 
is used to demonstrate the working procedure. Figure 2 is a photograph 
of the flight assembly and Table ©contains the input for the programs. 

The program provides hard copy computer prints of the projectile outline. 
Figure 3, and the corresponding estimated aerodynamic performance. Figures 
4 through 8.- Table 4 is a summary of the results. 

The available range data is superposed by an asterisk. 
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Figure 1. Outline of Long Rod Projectile 














































TABLE 1. SYNTHESIS OF DRAG COEFFICIENT EQUATIONS 
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TABLE 2. SYNTHESIS OF NORMAL FORCE AND STATIC MOMENT SLOPE COEFFICIENT EQUATION-' 
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TABLE 3. RETARDATION, ACCURACY FACTOR AND INITIAL YAW EQUATIONS 
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Accuracy Factor 

Ct I 

T _ La y 

i ,2 

Ma md 
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Initial Yaw Period 

Peak to Peak Distance 

As = IT 
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Figure 2. Photograph of 25/12 Projectile 
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Figure 3. Outline of 25/12 Projectile 

























TABLE 4. SUMMARY OF CALCULATED AERODYNAMIC COEFFICIENTS 
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Figure 4. Normal Force Slope Coefficient vs Mach Number 
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Figure 5. Static Moment Slope Coefficient vs Mach Number 
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Figure 6. Zero Yaw Drag Coefficient vs Mach Number 
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LIST OF SYMBOLS 


A 

b 

c 

c 

r 

c.g. 

c.p, 

d 

e 

g 

h/2 

j 

k 

£ 


£ 


grl 

£ - 
gr2 

A£ 


gr 


£ 


n 


o,a. 


m 

n 

s 

t 

V 


3 TAN 0 ) , operational parameter 
Intercept of vs M characteristic 
Slope of Cj^ vs M characteristic 
Fin blade length at root 
Fin blade length at tip 

Center of gravity of projectile, nose datum 

Center of pressure of normal force 

1.0 cal , reference diameter 

Base of Natural log 

Fin dimension (root recess) 

Fin dimension (blade height) 

Fin dimension (blade slant height) 

Fin dimension (blade extension from body) 
Cylindrical body length 
Groove length from nose (starting groove) 
Groove length from nose (last groove) 

Groove length 
Nose length 

Overall length of projectile 

Z + Z 
a n 

Mass of projectile 
Number of fin blades 
Range 

Fin dimension (average thickness) 

Velocity of projectile 
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Av 

X 

AR 


LIST OF SYMBOLS (continued) 

Velocity decrement over specified range 

Distance along projectile, nose datum 
2 

h /S„, Aspect ratio of fin planform 


A^ase body Area of body exposed to base pressure (cal)' 


base fin 


ref 


wetted body 


wetted fin 


Area of fin exposed to base pressure (cal)' 

2 

Reference area (.785 cal ) 

Area of body producing viscous flow drag 
Area of fin producing viscous flow drag 

Drag Force ^ zero-yaw drag coefficient 
P ^ \ef 


Sbb 

Sbf 

^DGR 

^DT 

^DTB 

‘"DTF 

^DVB 

^DVF 

^DWB 

^DWF 

^La 

Cw 

Ma 


Pressure drag coefficient - base of body 
Pressure drag coefficient - base of fins 
Drag coefficient due to grooves 
Total drag coefficient 
Total body drag coefficient 


Total fin drag coefficient 
Viscous drag coefficient - body 
Viscous drag coefficient - fins 
Wave drag coefficient - body (nose) 
Wave drag coefficient - fin 
Lift Force 


h P V 6 

Static Moment 
^2 P A^^^ d6 


, aerodynamic lift slope coefficient,6 = sin a_ 


, aerodynamic moment slope coefficient 
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LIST OF SYMBOLS (continued) 


^MaB 

^MaF 

Cx, 

Na 


Static moment coefficient - body 
Static moment coefficient - fins 


Normal Force 
>S P « 


aerodynamic normal force slope coefficient 



K 

M 

M, 

M 

Q 

R 


Normal force coefficient - body 
Normal Force coefficient - fins 
Axial moment of inertia 
Transverse moment of inertia 


I C, 
y La 

md^ ^Ma 


, aerodynamic jump factor 


Interference factor 

Mach number 

Mach number at muzzle 

Mach number along trajectory 

A » b 
D ret 

Operational parameter, ~ 2 m - 


c M + b 

Operational parameter, —- 

o 


a 

3 

6 

? 

A 

IT 


Angle of attack, employed here as a subscript 
2 1/2 

(M - 1) , operational parameter 

sine of total angle of attack 

Dispersion parameter, employed here as a subscript 
c^/c^, fin tip ratio 

Z.Ultb 
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LIST OF SYMBOLS (continued) 


P 


Ambient air density 


Fin sweep angle 
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APPENDIX A 


INITIALIZATION INSTRUCTIONS 

Initialization requires only the entry of the geometry and physical 
properties of the projectile, the range excursion and the muzzle Mach 
number as directed by the program listing according to the following 
schedule. 


Tektronix 

Figure 1 

Notes 

Nomenclature 

Identification 


L 

1 

n 

Refers to equivalent^^ conical 
nose length 

Li 

a 

a 

Includes cylindrical length of 
body to which fins are attached. 

A taper hub is approximated as 
a cylinder of half the length 
of the taper. A shallow taper 
aft body (BM6) requires an 
engineering evaluation. 

PI 

groove 

Distance to start of groove. Nose 
datum. 

P2 

AS, groove 

Length of groove set. 

H2 

h/2 

Radial fin length above hub. 

L 

G 

c 

r 

g 


K 

k 


T 

t 

Taken as maximum thickness of fin 
blade. 


Cl 

N 


n 


C2 

MO 

S 

W 

I 

I 1 
D 


c.g. 


M. 



I 


I 

d 


X 

y 


Restricted to either four or six 
blades in static moment calculation. 
Unrestricted in drag calculation. 
Nose datum. 

Muzzle.Mach number. 

Maximum range in meters 


May be given in either inches or 
millimeters. 
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FLOWCHART FOR AUTOMATIC PLOTTING ROUTINES 
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IDENTIFICATION OF FLOWCHART STEPS 


(1) Program Inputs: (All values are to be inputted in calibers unless 

otherwise noted) 

A - Number of Fins (4 or 6) 

B - Conical Nose Length 
C - Cylindrical Body Length 
D - Groove Location/Nose 
E - Groove Length 
F - Fin Dimension (Blade Height) 

G - Fin Blade Length at Root 
H - Fin Dimension (Root Recess) 

I - Fin Dimension (Blade Extention from Body) 

J - Fin Dimension (Maximum Thickness) 

K - Fin Blade Length at Tip 
L - Center of Gravity of Projectile 

M - 1.0 Caliber Reference Diameter (Millimeters or Inches) 

N - Normalized Weight of Projectile (Cal^) 

0 - Axial Moment of Inertia (Cal^) 

P - Transverse Moment of Inertia (Cal^) 

Q - Mach Number at Muzzle 
R - Maximum Range (Meters, < 4000) 

S - Plotting Device Number (l=Pen Plotter, 32=Screen) ■ 

(2) Projectile Design: Rear and Profile Views of Projectile and Fin 

(3) Optional Printout of Original Data 

(4) Program Selection: (Use User Definable Key) 

A - Nomenclature — A Listing of Aerodynamic Coefficients and their 

Definitions (Key 1) 

B - Table of Coefficients — A Listing of the Aerodynamic Coefficients 

and Definitions (Key 2) 

C - Total Normal Force Coefficient (C_ ) versus Mach Number Plot (Key 3) 

^ot 

D - Total Drag Coefficient (Cq^) versus Mach Number Plot (Key 4) 

E - Velocity versus Range Plot (Key 5) 

F - Total Pitching Moment about the Center of Gravity (C^^ , CG) versus 

Mach Number Plot (Key 6) 

G - Aeordynamic Jump Factor (J ) versus Mach Number Plot (Key 7) 

H - Reprint of Projectile Design (Key 8) 

I - Return to Menu (Key 9) 

(5) Stop 
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LONG ROD FINNED PROJECTILE DESIGN 
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ALL UALUES ARE IN CALIBERS UNLESS OTHERHISE NOTED 
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ENTER YES FOR HENU 
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PLEASE PRESS USER DEFINABLE KEY TO CONTINUE 










zn 

d t. 

O if 



31 

Xf M ** 



■d 

p p 



O 

if*4 if 



Xk 

— u. u 




•«4 



if 

4^ if if 




w»— jC 



•«4 

if .1^ •V 



** 

•->44 



u 

O U 44 



if 

1. if d 



•”> 

a—> o 



O 

o ja 



1. 

if 1. d 

ai 3 

¥1 

a 

z a 

d d 



44 44 44 

O O 

31 ¥k 

if 

if P P 

St 91.0 

-O m WH. p 

JZ 3f !.£ if if 

d 

O C ^ 

4-»7S *<4 

if O if 

11^4* 


0'4> P U if U 

— JQ — 


i.ja 1. if->4 c-^ 



o 

4-^ 

44 «f44 

f-. 11 114^ 


If f:'«- o'*- W'*- 
— Of if ¥1 «l 
«-• •V M> o if O 
e*# u PM-o 1.0 
if u-^ if ii a 


u*— U u--« 

if.^ QfCf>*^-r4-^ if+# 

*"> O u U 

O U if U U O if 

C if 1. if if 1. ■ 


LU 

•<.4 if -fH O 44 44 

a-> O'*- o •■*>—> CL o 

0 ^ 

U'-dH. uO P'4- P 

O O 1. O O 1. 

3 

o if <4 if o if 

ifcif ai.i.ifa 

H- 

1 - O'*- if c r. 

£ ojc p a ax; 

o: 

'*- ao'*- u o 1 . o 

44 44 O if 44 if 


if d 1 .Z ifz 

if •> 4 x: if if z 

o 

O il 44 O O 44 

C£ 1.44 44 jC£ 1.44 

Z ' 

Oi= poll. 91 p m 

O-*^ O u 4444 0 

Ui 

44 If P if P 

N*. '*- if 1. '*- 1. 

z 

if C if^-tH U *4 

1 . M 0 1 . 1 . 0 

o 

u 1 . o u dJ3 z 

en o**' '*- 0 o-*4'*- 


1. OZ L. C u Irt u 

d'*- Pd '*-'*- P 


O'*- O 1.44 344 

C if if 44 if 44 


tj. 9)U. O’lH (Z-mH 
** P XCL-^O. 
*— P*^— c 
0 VjC d-"- — 

u c o d 7) d 

i, U4^ 

O *'4 o o o o 

Z'*‘Q.ZI->- >H* 

d 

L. if 

9>jC 


if if if if if if 
£ OXlXiJCJZ 
•M u 


PuOifPPuif 
V if-^ o -^ if 

><tH*4> 4. u 

d u'+’ u uH-.tH 

3*^ if ifs*. 

u o ifo if 

o if o if if o 

09)9)u009tu 

u d d u u d 


4-'7f p 1.1. m 


4. QfH.4>4>H. OH- 

O 1 . o o o o o 
d c. 

ifififififififif 

aiA a Q. a o.** cl 

OifOOOOPO 

(n Q. (o (f) in o) o (f) 

L3 

o 

CL 

COflQCQli-H‘f~Of“ 

<ra.<c-cco:<x f <c 

ZOZZZZL9Z 

0X000000 


w 

if Wd ”d d W Wd d 
•>^ ••H d 1. d d 1. 
> u 1. Ii ifd 1. C Wd 
d-4-od*- d*dd 

o-4>— O — 

m'4. Qfu^dififud 
iiiniio-»-'>irt«rt4^ 
H-od*<^i.Odd-'4 0 
OOCD 3 >CLf-Za>S-h- 


a)OQma:cQu.u.u.iju 

3 CO 3 O f-3 OQ 3 h- 

ooooooooo 


39 


CDT Total drag coefficient 

CLA Slope of the Lift Coefficient 

J ZETA Aerodynanic juHp factor 


STATIC AERODYHAHIC COEFFICIENTS FOR LONG ROD FINNED PROJECTILES 
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HACH No 





CDT vs. MACH No 
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UELOCITY vs. RANGE 
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CMATyCG vs. MACH No 
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J ZETA vs. MACH Ho 
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